Manifestations of magnetic vortices in the equation of state of a Yang-Mills plasma 
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The vacuum of Yang-Mills theory contains singular stringlike objects identified with center (mag- 
netic) vortices. Percolation of magnetic vortices is known to be responsible for the color confinement 
in the low-temperature phase of the theory. In our work we study properties of the vortices at finite 
temperature using lattice simulations of 5(7(2) gauge theory. We show that magnetic vortices pro- 
vide a numerically large contribution to thermodynamic quantities of the gluon plasma in Yang-Mills 
theory. In particular, we observe that in the deconfinement phase at temperatures T c < T < 3T C 
the magnetic component of the gluon plasma produces a negative (ghostlike) contribution to the 
anomaly of the energy-momentum tensor. In the confinement phase the vortex contribution is 
positive. The thermodynamical significance of the magnetic objects allows us to suggest that the 
quark-gluon plasma may contain a developed network of magnetic flux tubes. The existence of the 
vortex network may lead to observable effects in the quark-gluon plasma because the chromomag- 
netic field of the vortices should scatter and drag quarks. 

PACS numbers: 12.38.Aw, 25.75.Nq, 11.15.Tk 



I. INTRODUCTION 

Studies of properties of thermal plasma became a ma- 
jor development in QCD in recent years, for a review 
see, e.g. PJ,[2J]. Properties of the plasma are studied both 
directly, at RHIC and via lattice simulations. On the 
theoretical side, novel ideas, like AdS/CFT correspon- 
dence are being invoked [3| , to say nothing of traditional 
approaches based on various quasiparticle models [J] and 
on field theory at finite temperature. 

The traditional approach to the thermal plasma treats 
it, in zero approximation, as gas of free gluons and quarks 
and, then, takes into account perturbative corrections. 
An outcome of such calculations is a representation of 
the energy and pressure densities as perturbative series 
in the effective coupling constant g 2 (T): 



e(T) = CsB T 4 /pert(5 2 (T)), 



(1) 



where csb is the Stefan-Boltzmann (SB) coefficient, pro- 
portional to the number of degrees of freedom, and the 
perturbative expansion, / pe rt(<? 2 (T)) is known explicitly 
up to terms of order g 6 lnT. Perturbative predictions for 
such global characteristics of the plasma as energy den- 
sity turn to be in reasonable agreement with the data [2| • 

On the other hand, some particular properties of the 
plasma, such as viscosity [fj|, indicate that, in the zero 
approximation, plasma is to be considered rather as an 
ideal liquid than an ideal gas. There is no yet a coherent 
picture that would unify both perturbative and nonper- 
turbative features of the QCD plasma. 

It was speculated in Refs. |g, [D, [8( that there exists a 
magnetic component of the Yang-Mills plasma which at 
temperatures not much above the critical temperature 



T c is crucial for the plasma properties. In Refs. 0, Q 
constituents of the magnetic component arc thought to 
be classical magnetic monopoles. In Ref. Q the magnetic 
component is identified with so-called magnetic strings 
related to magnetic monopoles. The properties of the 
strings, or center vortices and their role in confinement 
have been discussed in the lattice community for more 
than a decade, for review and references see Ref. [S|. 

According to the vortex picture the quark confinement 
emerges due to spatial percolation of the magnetic vor- 
tex strings which lead to certain amount of disorder. The 
value of the Wilson loop changes by a center element of 
the gauge group if the magnetic vortex pierces the loop. 
Therefore, very large loops receive fluctuating contribu- 
tions from the vortex ensembles. These fluctuations make 
the average value of the Wilson loop very small. One can 
show that the suppression of the loop follows an area law 
for very large loops |i|. 

Magnetic monopoles are certain gluonic objects that 
are related to color confinement via the so-called dual 
superconductor mechanism [lOJ . The condensate of the 
monopoles - which is formed in the low-temperature 
phase - expels the chromoelectric field exhibiting a dual 
analogue of the Meissner effect. The chromoelectric 
field of quarks is squeezed into a dual analogue of the 
Abrikosov vortex leading to quark confinement (for a re- 
view see Ref. |11|). 

The presence of both the magnetic monopoles and the 
magnetic vortices in the vacuum of Yang-Mills theory 
seem to be an indication of the existence of a more com- 
plicated genuine non-Abelian object. Indeed, Abelian 
monopoles and center vortices appear to be strongly cor- 
related with each other [12]: almost all monopoles are 



sitting on top of vortices. In SU(2) gauge theory the 
genuine object is a monopole- vortex chain 12]. In SU(3) 
Yang-Mills theory one can expect existence of monopole- 
vortex 3-nets: the junction of the three vortices ma y b e 
viewed as a nexus [13[ or/and as a center monopole [lj] . 
An illustration of the chains and nets in a 3-dimensional 
timeslice is shown in Figure [T] 




SU(2) chains 



SU(3) nets 



plasma. Preliminary results of our investigation were re- 
ported in Ref. [231 . Our lattice calculations refer to the 
case of pure SU{2) gauge theory which is conceptually 
similar to a more realistic theory with three colors. 

The structure of the paper is as follows. In Section HT1 
we discuss basic relations of Yang-Mills thermodynamics 
both in the continuum limit and on the lattice. In Sec- 
tion IIIII we present our results for the thermodynamics 
of SU(2) Yang-Mills theory on moderately-sized lattice. 
Although the numerical results on the SU(2) thermody- 
namics are available in the literature 24\ we repeat the 
calculations on lattices of particular geometries as these 
lattices are later used in Section [IV] to calculate ther- 
modynamical properties of magnetic vortices such as the 
thermal vortex density and contribution to the (trace) 
anomaly of the energy momentum tensor. A short sum- 
mary of our results is given in Sect ion [V] which is followed 
by a concluding Section. 



FIG. 1: (Color online) An illustration of the monopole-vortex 
chains and the monopole-vortex nets in SU(2) and SU(3) 
Yang-Mills theories, respectively. 

The formation of these composite objects is essential 
for the self-consistent treatment of monopoles in the 
quark-gluon plasma [6| . Similar monopole-vortex chains 
were found in numerous (non-)supersymmetric non- 
Abclian gauge theories involving various Higgs fields [la, 
llq ]. In particular, in Abelian gauge models with a com- 
pact gauge field, monopoles and vortices should be orga- 
nized in chains/nets in order to be consistent with break- 
ing of the confining electric string in specific representa- 
tions of the gauge group [lg|. A short overview of the 
monopoles and vortices in the Yang-Mills plasma can be 
found in Ref. p7| . 

In terms of the continuum theory, the magnetic strings 
are defined as (infinitely thin) surfaces that are closed in 
the vacuum state and can be open on an external 't Hooft 
line (the same as the standard, or "electric" strings can 
be open on the Wilson line). One can argue (for review 
and references see Ref. |18|), that the magnetic strings 
of the continuum theory can be identified with the cen- 
ter vortices of lattice studies. Moreover, it is also known 
that, indeed, properties of the vortices are changing dras- 
tically once temperature is increased above the critical 
value T c H, [H [13, [H H3]. In particular, they become 
time oriented and this agrees with the assumption that 
the vortices become a component of the Yang-Mills ther- 
mal plasma [6J. 

Once the magnetic component of the plasma is iden- 
tified with the magnetic vortices, further information on 
its properties can be obtained by direct numerical cal- 
culations on the lattice at finite temperature. In this 
paper, we report on results of lattice calculations of con- 
tribution of magnetic strings to the trace anomaly of 
energy-momentum tensor in the gluon plasma. The trace 
anomaly is intimately related to the equation of state, or 
energy density and pressure of the thermal Yang-Mills 



II. YANG-MILLS THERMODYNAMICS 
AND TRACE ANOMALY 

A. Continuum thermodynamics 

The free energy of an SU(N) gauge system 

F = -T log Z(T,V), (2) 

is expressed via the partition function Z as follows: 



Z = |DAcxpj-^TrG^j 



(3) 



where G^ v = G a „ u t a is the field strength tensor of the 
non-Abelian field A, and t a are the SU(N) generators 
normalized in the standard way, Tr t a t b = \5 ah . The 
pressure p is given by the derivative of the partition func- 
tion © with respect to the spatial volume, V, 



P 



T dlogZ(T,V) _ F _T 
V d logF " ~V ~ V 



log Z(T,V). (4) 



The last two equalities are valid for a sufficiently large 
and homogeneous system residing in thermodynamical 
equilibrium. 

The energy density e is given by the derivative of the 
partition function ([3]) with respect to the temperature T, 



e = 



T d log Z(T,V) 
V dlogT 



(5) 



The relation between the pressure (HJ and the energy ([5]) 
constitutes the equation of state of the system, p = p{e). 
According to Eqs. ^ and ([5]) it is sufficient to de- 
termine the partition function of the system in order to 
calculate the equation of state. However, available lat- 
tice simulations techniques are suitable for the calcula- 
tion of quantum averages of operators rather than for the 



evaluation of the partition function itself. On the other 
hand, both the energy density e and the pressure p can 
be derived from the quantum average of a single quantity, 
which is the trace of the energy-momentum tensor T^. 
In SU(N) gauge theory the energy-momentum tensor 
is given by the formula 



T„ 



2Tr 



1 



^ fxcr^va a ^fiu^ap^ap 



(6) 



At the classical level the energy-momentum ^ is trace- 
less because the bare Yang-Mills theory is a confor- 
mal theory. However, at the quantum level the confor- 
mal invariance is broken, and, consequently, the energy- 
momentum tensor exhibits a trace anomaly: the quan- 
tum average of the trace of the energy-momentum tensor, 



0(T) = (Tjf) 



3p, 



(7) 



is nonzero. The trace anomaly is intimately related to 
the gluon condensate which breaks the scale invariance 
of the theory 



where the Gell-Mann-Low /3-function is 
Of, 



0(9) 



-g 6 (b + hg z + ...). 



dlog^ 
with one- and two-loop perturbative coefficients 



bo = 



UN 
3(4tt) 2 



h = 



347V 2 

3(4tt)' 



(8) 



(9) 



(10) 



The thermodynamic relations in Eqs. {4| and |5|) relate 
the trace anomaly ([7]) with the pressure p, 

Then, the pressure and, respectively, the energy density 
can be expressed via the trace anomaly, as follows: 



a f T dTi 6>(Ti) 

s(T) 



3T q 



dT\ g(Ti) 



0(T). 



(12) 
(13) 



Equations (fT2"|) and (fT3")) demonstrate that the trace 
anomaly ([7]) is the key quantity that allows us to re- 
construct the whole equation of state. 

The trace anomaly should vanish in a system of free 
massless relativistic particles 



-SB 



3p 



SB 



N d f — T 4 
30 



(14) 



[where iVa.o.f. = 2(7V 2 — 1) is the number of degrees of 
freedom in the noninteracting gas of the SU(N) gluons] 



or in the case when excitations are too massive com- 
pared with the temperature, m 3> T (then e ~ p ~ 
exp{— m/T}). For Yang-Mills theory these statements 
imply that the dimensionless quantity 9 /T 4 should ap- 
proach zero at very high temperatures (the gluons form 
a weakly interacting gas) and at very low temperatures 
(note that in SU(N) Yang-Mills theories the mass gap 
is much greater than the critical temperature }27l|). The 
latter property is used to fix the lower integration limits 
in Eqs. ||Ij2J and {B} at T = 0. 



B. Lattice thermodynamics 

The lattice analogue of the partition function ^ of 
SU(N) gauge theory is 



Z(T,V)= /cC/exp{-/3^5p[C/]} 



(15) 



The plaquette action Sp [U] of the gluonic link fields U xll 
is usually written in the Wilson form, 



SplU] = 1-— ReTrC/p. 

L J N 



(16) 



The identification of the lattice results with the physics 
in continuum is achieved in the limit of vanishing lattice 
spacing, a — > 0: 

x+aft 

U XI1 = expjig / dyA^y)} (17) 

X 

-> 11 + iagA^x) + 0(a 2 ) , 
The lattice spacing a is a function of the lattice coupling 



& = 2N/g 2 



(18) 



The spatial volume V = {L s a) 3 and the temperature 



1 
L t a 



(19) 



of the gluonic system are related to the asymmetric geom- 
etry of the Euclidean lattice, L? s L t . The shorter direction, 
L t with L t ^ L s , is associated with the imaginary time, 
or, "temperature" direction. The imaginary time formal- 
ism allows us to calculate various thermodynamic quan- 
tities corresponding to the gauge system residing in ther- 
modynamic equilibrium at given temperature and vol- 
ume. Using the relation T(d/dT) = —a(d/da) one can 
rewrite the continuum expression for the anomaly (|11[) in 
a form suitable for numerical lattice simulations, 



gen 

rpA. 



= 6L1 



d /3(a) 
d log a 



({Sp)t - (Sp)o) ■ (20) 



Here the plaquette averages (Sp) T and {Sp) Q are the ac- 
tion densities taken, respectively, in a thermal bath at 



T > (L 3 s L t lattices) and at T = (L 4 lattices). In 
Eq. (J2TJ)) it is implied that the T = plaquette expec- 
tation value is subtracted to remove the effect of quan- 
tum fluctuations, which lead to an ultraviolet (UV) di- 
vergency of the quantum expectation value. As a result, 
the trace anomaly becomes a UV-finite quantity, which 
is normalized to zero at T = because of the existence 
of the mass gap. Equation (|2"0)) is a lattice version of 
Eq. ([8]). The trace anomalies and the equation of state 
for SU(2) and SU(3) gauge theories were calculated in 
Refs. |24| and [25|], respectively. 

The anomaly P0|) can be separated into electric and 
magnetic parts, respectively, 



e(T) = e E (T) + e M (T) 



(21) 



The electric contribution comes from the field strength 
tensors of the chromoelectric fields, Ef = Gf 4 , while 
the magnetic part is solely due to chromomagnetic fields, 
Bf = (l/2)ejkG^ k . On the lattice the former is associ- 
ated with temporal plaquettes, P t = Pa with i = 1, 2, 3, 
while the latter is related to the spatial plaquettes, P s = 
Pii with i<j = 1,2,3: 



1E 



(T) 



y4 

9m(T) 



3L 4l'dP(a 



3L 



• l a log 

i i 9P(a 



a log 



i) ■ (<&,>, 



(Sp)a) , (22) 
(S P )o) ■ (23) 



Summarizing, the only lattice observable which is 
needed for the calculation of the trace anomaly (f^TJj) [in- 
cluding its electric (|2"2"j) and magnetic (f2"2")l parts] and, 
consequently, for the determination of thermodynamics 
of Yang-Mills theory, is the difference between the expec- 
tation values of the action densities calculated at zero- 
temperature (at L t = L s ) and at finite temperature (at 
L t < L s ). The /3-function in Eq. (j20|) has to be deter- 
mined from independent lattice simulations. 



III. SU(2) THERMODYNAMICS AT 
MODERATE LATTICES 



In this Section we present basic results for the ther- 
modynamics of SU(2) gauge fields. The thermodynam- 
ics of SU(2) lattice Yang-Mills theory has been studied 
in detail in Ref. [2J]. In this Section we repeat certain 
numerical calculations of Ref. [2J] at moderate sized lat- 
tices which will be used in the subsequent Section in the 
investigation of the vortex contributions to the thermo- 
dynamics of the system. 

We are working on the lattices 18 3 x 4 and 18 4 cor- 
responding to finite- and zero-temperature cases, respec- 
tively. At the asymmetric lattice with L t — 4 the thermal 
phase transition is realized at j3 = j3 c ~ 2.3, Ref. [20 |. 
In physical units the transition in SU{2) gauge theory 
is achieved at temperature [23] T c = 305(8) MeV, if one 
fixes the scale by setting the T = tension of the chromo- 
electric string at the phenomenologically accepted value 
of a 1 ' 2 = 440 MeV. 
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2.27 


2.3 


2.35 


2.4 


2.45 


2.53 


2.66 


T/T c 


0.909 


1.005 


1.188 


1.401 


1.651 


2.14 


3.24 


a, fm 


0.183 


0.165 


0.140 


0.119 


0.101 


0.0777 


0.0513 


L s a, fm 


3.293 


2.977 


2.520 


2.136 


1.813 


1.399 


0.9234 


d/3 
d log a 


0.2965 


0.2982 


0.3011 


0.3038 


0.3064 


0.3103 


0.3160 



TABLE I: Parameters of our numerical simulations: temper- 
ature (in units of the critical temperature T c ), lattice spac- 
ing a, and spatial extension L 3 a for selected values of the 
lattice coupling j3. 



The gauge field configurations are generated using the 
standard Wilson action (|16p . In order to determine the 
lattice spacing a as a function of the lattice gauge cou- 
pling /3 we use the interpolating function of Ref. [28| 
which describes with a high accuracy various lattice cal- 
culations of a. The lattice spacing a in units of the T = 
string tension a is parameterized as follows: 



2 647T 4 6o h 647r 4 b 

ln(aa) = — — + ^ ln ^ 



3&o/3 



(24) 



The first two terms in the right hand side (r.h.s.) of the 
above equations represent, respectively, the first and the 
second order loop corrections coming from the perturba- 
tion theory. The coefficients bo and b\ are the coefficients 
of the /3-function ^ which are given by Eq. pH]) with 
N = 2. The other two terms in the r.h.s. of Eq. (|2~4")) are 
introduced to mimic higher-order perturbative correc- 
tions as well as nonperturbative effects. It was found [28| 
that the choice c = 4.38(9) and d = 1.66(4) reproduces 
measured valued of the lattice spacing a very well. 

We use the interpolating relation (|24[) in order to de- 
termine the temperature (| 19[) and the prefactor in the 
r.h.s. of Eq. (J2"D)) . In Table U we present several impor- 
tant quantities inherent to our numerical calculations. 

In Fig. [2] we show the trace anomaly PD]) along with 
its electric (|22|) and magnetic (|23|) parts. It is very in- 
teresting to notice that the magnetic contribution to the 
anomaly vanishes in the deconfinement phase at a tem- 
perature approximately twice larger than the transition 
temperature, 



9m(T )=0, To = 2.15(1) T c 



(25) 



This distinguished value of the temperature is also seen 
from the asymmetry of the so-called A 2 -condensate: the 
thermal fluctuations of the electric and magnetic com- 
ponents of the gluonic fields in the Landau gauge be- 
comes equal around 2T C , see Ref. [29j for details. In the 
case of the ST/ (3) gauge group the magnetic contribution 
to the trace anomaly also vanishes around the temper- 
ature which is slightly higher than 2T C , Ref. [25]. The 
mentioned coincidence allows us to conjecture that the 
value (|2"5")) of the critical temperature, T rs 2T c , is in fact 
universal for all numbers of colors N. 




T/T 



FIG. 2: (Color online) The trace anomaly 9, represented as a 
dimensionless ratio 6/T 4 , as a function of the temperature T 
(in units of the critical temperature, T c ). We show the full 
anomaly (circles), and its electric (squares) and magnetic (tri- 
angles) contributions determined by Eqs. (f20|) . (f22|) and ([23l) . 
respectively. 



IV. THERMODYNAMICS OF VORTICES 



Magnetic vortices and gluon trace anomaly: 
lattice definitions 



As we have discussed in the Introduction, (magnetic) 
center vortices are stringlike configurations of the gauge 
fields which are relevant for color confinement in the low- 
temperature confinement phase. In this Section we inves- 
tigate how these objects contribute to the thermodynam- 
ics of the SU{2) gauge system. 

The magnetic vortices carry a magnetic flux character- 
ized by a nontrivial center element of the gauge group. 
Such objects can conveniently be found in a specific 
gauge, which makes the gauge field as close to the cen- 
ter subgroup of the gauge group as possible. In 517(2) 
Yang-Mills theory the relevant gauge - often called as 
the (direct) maximal center gauge - is defined by a max- 
imization of the functional 



maxFf[/ s 

Q 



F[U] 



i 



(TrUtf 



with respect to the gauge transformations 



U 



Xfl 



Xfi 



x xfx x-\-£l ' 



(26) 



(27) 



In Eq. (j26|) the sum goes over all links / of the lattice. 

The center element of the gauge field Ui can be identi- 
fied after the SU{2) maximal center gauge (|26|) is fixed: 



Z, = signTr[/, = ±l 



Zi ez 2 



(28) 



The maximal center condition (|26[) fixes the gauge free- 
dom (f27|) up to a center subgroup of the SU{2) gauge 
group: 



U X4 i —* U x — uj x U xfi uj x+ f l , 



uj x e Z 2 



(29) 



In terms of the Z 2 gauge variables (|28]) the center trans- 
formation ([29|) . 



Z 



Xfl 



^XfJ. — WxZxflUx + fL , 



(30) 



can naturally be interpreted as a Z 2 transformation in a 
Z 2 effective gauge theory written in terms of the gauge 
fields {HP. 

The field strength tensor of the center gauge fields ([28]) 
is the Z 2 plaquette 



Zp — Z\ZiZ^Z± 



Z P £l 2 , 



(31) 



where the subscripts 1 ... 4 denote the links forming the 
border of the plaquette P. The Z 2 -field strength ten- 
sor pip is invariant under the Z 2 gauge transforma- 
tions (j30l) . 

The vortex worldsheets S are identified with the help 
of the Z 2 -plaquette (f3Tj) . Let us introduce the notation 
*P for the plaquette which belongs to the dual lattice and 
which is dual to the plaquette P. Then the dual plaquette 
*P does not contain the center vortex if the correspond- 
ing plaquette is center-trivial, Zp = +1. Equivalently, 
one may say that no vortex is going through the plaque- 
tte *P, if Zp = +1. However, if Zp = —1, the plaquette 
*P is a part of the vortex world-sheet. In short, 



Z P = 



-l, if Pes 
-1, if P <£ S 



(32) 



One can easily prove that vortices are closed loops since 
their worldsheets are closed surfaces. 

Using Eq. (|32|) one can define the lattice vortex density 



Pp = ^i l - z p) 



1, if P e S 
0, if P 4 S 



Then the averaged lattice vortex density is 

P=\(1-(Z P )). 



(33) 



(34) 



The lattice fields can be decomposed into the vortex 
singular contribution and the rest 

Ui =Z t Ui, (35) 

where the 517(2) gauge field Ui is defined as follows 

Tr Ui > . (36) 

The trace of the SU(2) plaquette (which is an explic- 
itly 5(7(2) invariant quantity), can be decomposed into a 
singular (center-valued) contribution and a regular con- 
tribution, respectively, 

Tr Up = Z P Tr U P , U P = Ux U 2 U% U\ . (37) 

This relation is very useful for us since it allows us to 
identify the contribution of the singular magnetic vor- 
tices into the vacuum expectation value (vev) of the 



SU(2) plaquette (Tr Up). This quantity is related to 
the vev of the gluonic action (|16[) . which, in turn, con- 
tributes to the thermodynamics of the system via the 
trace anomaly (f2"0f . 

In order to calculate the contribution of magnetic vor- 
tices into the trace anomaly (|20[) . and, consequently, into 
the equation of state, we utilize the following chain of 
considerations. For any given configuration of the gauge 
field Ui one can identify the location of the center vortex 
worldsheet, *£, as we described above. The worldsheet 
*£ is a collection of plaquettes *P of the dual lattice, 
which forms a closed surface, FigureO Then the set of all 




FIG. 3: (Color online) A three dimensional visualization of 
the magnetic vortex string * E which pierces a set of the pla- 
quettes PgE. 

plaquettes of the lattice can be separated into two (gen- 
erally, unequal) subsets: (i) the plaquettes P of the orig- 
inal lattice belonging to the vortex worldsheet, P G S; 
and (ii) the plaquettes P which are outside of vortices, 
P $_ S. The total action of the gauge field configura- 
tion can equivalcntly be decomposed into the sum of the 
contributions of the individual plaquettes belonging to 
the vortex worldsheet and the sum going over plaquettes 
outside of vortices, respectively: 

J2s p = J2 s p + J2 s p ( 38 ) 

p Pes p^s 

= ^p P Sp + ^(l-pp)Sp. (39) 
p p 

Then we divide both sides of Eq. (|3"5|) by the total num- 
ber of the lattice plaquettes, 6L^L t , and average over all 
gauge field configurations. As a result, we get that the 
average gluonic action per plaquette is given by the sum 



(S P ) = (S P ) von + (s P y 



(40) 



where (Sp) VOIt is the vev of the gauge action coming from 
the vortex worldsheets, 

(Sp) vmt = ( Pp Sp) = -^ i -(J2ppSp) (41) 



p+-((TrU P )-(TrU P ) 



and the action coming from the space unoccupied by the 



vortex worldsheets, 

{S P ) aovmt = ((I-pp)Sp) 

* <£(l-pp)5 P ) 



(42) 



mu t 



i 



l-p--[{TrUp) + {TrU P ) 



In order to derive the relations Eq. (|41|1 and (T42j) we 
used Eqs. (J33J) , flM} as well as Eq. HH) for the case of 
two colors, N = 2. Then we utilized the relations, 

p P Sp = p P + \z p Tt Up - l -Tr Up , (43) 

(I-Pp)Sp = (1 - p P ) - ^Z pTr Up -^Tr Up, (44) 

and used the fact that, according to Eq. (|3~7| . 



Z P Tr Up = Tr U P . 



(45) 



It is important to realize that the vev of the quan- 
tity (|45p can be interpreted as the vev of the plaquette 
Tr Up evaluated on configurations with "removed vor- 
tices" . Indeed, if we manually set Zp ~ 1 at all plaque- 
ttes of the gauge field configuration, then the vacuum 
expectation of Tr Up at these modified configurations is 
automatically equivalent to the vev of Tr Up at the orig- 
inal (unmodified) configuration: 



(Tr Up) = (ZpTrUp) = (TrU P ) T 



emoved vortices 



(46) 



The procedure of "removing vortices" is equivalent to a 
shift of the gluonic link fields Ui , 



Ui -> Zi Ui = Ui , Z t =signTr J7, 



(47) 



according to the decomposition (|35|) . By definition, 
the thus-modified configuration contains no center vor- 
tices. Note that the gauge field configurations with re- 
moved vortices are not able to support the color confine- 
ment 30]. 

Summarizing, the lattice density of the gluonic ac- 
tion - which is crucial for the determination of the trace 
anomaly (|20|) - can be split according to Eq. ([40|) into a 
piece coming from the vortex worldsheets (|41j) and arising 
from the regions of space lying outside of vortices (|42[) . 
In the next Subsections we analyze these contributions 
in more details. 



B. Center vortices and action density: 
numerical results 



We have studied vortex properties on the lattices 18 4 
and 18 3 x 4. Depending on the value of the lattice cou- 
pling (3 we used from 100 to 500 statistically independent 
configurations for the zero temperature lattice 18 4 , and 
from 700 to 1600 configuration for the finite-temperature 



lattice 18 3 x 4. Generally, the larger /3 the larger is the 
number of gauge field configurations needed to reach an 
acceptable statistical accuracy of the numerical results. 
The maximal center gauge (|26l) was fixed using a sim- 
ple iterative algorithm. In our simulations we considered 
one Gribov copy only. Despite the unsophisticated na- 
ture of our approach, we believe - following the results 
reported in Ref. 31] and in Ref. [32| - that this algorithm 
allows us to capture not only qualitative but also, to a 
large extent, various quantitative features of the vortex 
ensembles at finite temperature. In our visualization of 
the vortex-related data we indicate the statistical errors 
only, leaving aside the (systematical) error bars related 
to the gauge fixing issues such as the choice of the Gribov 
copy. 

Below we present results of our numerical simulations 
for several important quantities calculated both in lattice 
and in physical units. Before proceeding to the calcula- 
tion of the trace anomaly and other related quantities, we 
discuss the important numerical results for the density of 
the magnetic (center) vortices. 



1. Vortex density 

At finite tem per ature the vortex density p was studied 
in Ref. [H, [H HH, HI- We calculate the vortex density 
both at finite temperature and at zero-temperature in 
order to disentangle the finite-temperature effects from 
the effects related to the variation of the UV cutoff. In- 
deed, the UV cutoff - the role of which is played by the 
(inverse) lattice spacing a - is controlled by the lattice 
coupling constant (3, as, for example, illustrated by the 
interpolating function (|24|) . The temperature is related 
to the extension of the lattice in the short direction (fi"9|) . 
The comparison of the lattice results obtained at differ- 
ent lattice geometries and at the same fixed value of (3 
allows us to single out effects of temperature variations 
on the vortex density. 

In Figure 0] we show the total density of the center 
vortices (15"4")) calculated on the lattices 18 4 and 18 3 x 4. 
One can observe that in the confinement phase, (3 < (3 C ~ 
2.3, the vortex density is almost the same for L t = 18 and 
L t — 4 lattices. 

In the deconfinement phase, (3 > (3 C on the L t = 4 
lattice, the vortex density starts to deviate from its zero- 
temperature value (calculated on the L t — 18 lattice). 
It is quite surprising to notice that just above the phase 
transition the density of the thermal vortices gets lower 
compared with the value of the vortex density at zero- 
temperature. However, as the temperature (or, equiva- 
lently, the value of the lattice gauge coupling (3 on the 
L t = 4 lattice) increases further, the absolute value of the 
difference in the vortex densities at Lt = 4 and L t = 18 
lattices reaches its maximum, and then the deviation be- 
tween the vortex densities starts to decrease again. Thus, 
the effect of the temperature on the vortex density is 
quite nontrivial. 




FIG. 4: (Color online) The lattice density of the center vor- 
tices (|34l) at zero-temperature (the solid line) and at finite 
temperature (dotted line) vs the gauge coupling f3. The sim- 
ulations are done on 18 and 18 x 4 lattices, respectively. 
The vertical dashed line marks the critical coupling /3 C corre- 
sponding to the finite-temperature phase transition. On the 
18 3 x 4 lattice, a lower value of (3 corresponds to a lower tem- 
perature. Hereafter some sets of the numerical data points 
are connected by lines to guide the eye. 



At zero-temperature the vortex ensembles typically 
consist of one very long vortex worldsheet which spreads 
over almost the whole spacetime, and a lot of short-sized 
clusters of vortex worldsheets. The worldsheet of the in- 
finitely long vortex trajectory is usually associated with 
a condensed (infrared) vortex component while the short 
UV clusters form a perturbative (UV) component. It is 
the condensate of vortices which leads to the confinement 
of color charges. At T > T c the spatial string tension in 
the deconfinement phase is supported by the vortex per- 
colation in the spatial dimensions [Hj, [2l| . 

At finite temperature it looks natural to separate the 
strings in thermodynamic ensembles into two compo- 
nents describing, respectively, the virtual strings and the 
real strings. It is easy to explain the meaning of "vir- 
tual" and "real" in terms of particlelike objects which we 
understand much better than vortices which are string- 
like defects [6|: the real objects are thermodynamically 
relevant while the virtual ones are not. In other words, 
the virtual objects are associated with the ground state 
of the theory realized at zero-temperature while the real 
particles are the excitations of this ground state. 

Consider, for example, the Abelian monopoles at finite 
temperature in Euclidean Yang-Mills theory. At zero- 
temperature a part of the monopoles is condensed (simi- 
larly to vortices) while another part of the monopole tra- 
jectories belong to ultraviolet clusters. At finite temper- 
ature the ensembles of the monopole trajectories should 
contain both the virtual monopoles belonging to the vac- 
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uum and, at the same time, the thermodynamically rel- 
evant (or, real) monopoles. In order to measure physi- 
cally meaningful observables in the continuum limit, one 
should be able to separate these types of monopoles from 
each other. According to Ref. Q the thermodynamically 
relevant monopoles can be distinguished from the virtual 
monopoles by a simple principle: the thermal monopoles 
should wrap around the temperature (compactified) di- 
rection of the Euclidean space. Moreover, the quan- 
tum density of the thermal monopoles is not equal to 
the average density of the total length of the (wrapped) 
monopole trajectories as one could naively guess on gen- 
eral grounds. Such a quantity is divergent in the UV 
regime making its interpretation somewhat obscure. The 
density of the thermal monopoles corresponds to the den- 
sity of the winding number s of the monopole trajectories 
in the temperature direction of the Euclidean spacetime: 



thermal 
rmon 



1 



(N). 



(48) 



Coming back to the magnetic vortices, one can ask the 
important question: how to separate the thermal compo- 
nent of the vortex density from the density of the virtual 
component? In other words: in Euclidean lattice simula- 
tions at finite temperature we observe a set of closed vor- 
tex trajectories in each configuration of the gauge fields. 
Which vortex worldsheets are real (thermal) and which 
worldsheets corresponds just to virtual strings? Follow- 
ing our experience with particlelike monopoles [6| one can 
suggest that the worldsheets of the real vortices are char- 
acterized by a nontrivial wrapping number with respect 
to the (compactified) temperature dimensions. 

A vortex analogue of Eq. (|4"8|) remains to be derived. 
Note that the naive suggestion for the thermal vortex 
density - the averaged area of the wrapped vortex world- 
sheets per unit Euclidean volume - in general should not 
be correct. In fact, in Ref. [y] it was shown that the anal- 
ogous quantity for the monopoles (the average length of 
wrapped monopole trajectories per unit Euclidean vol- 
ume) is - contrary to the correct expression (|48|) - diver- 
gent in the UV limit and thus cannot serve as the density 
of the real monopoles. 

Let us stress again a somewhat puzzling behavior of the 
vortex density, Figure |U Naively, in a quantum field the- 
ory an expectation value of a local quantity at finite tem- 
perature should receive contributions from both quantum 
field fluctuations - which are inherent solely to the zero- 
temperature ground state - and from various thermal ex- 
citations. The (virtual) quantum field fluctuations have 
the same structure both at zero and at finite temperature 
(for example, the UV divergencies both at zero and at fi- 
nite temperatures are identical [331]). Therefore, in order 
to figure out a contribution of the thermal excitations 
one can fondly subtract a zero-temperature vev from the 
thermal expectation value. This subtraction, for exam- 
ple, is routinely performed for the trace anomaly of the 
Yang-Mills fields pp|) . The puzzle is that if the subtrac- 
tion procedure is applied to the vortex density, then the 



would-be thermal vortex density takes a negative value in 
the deconfinement phase, Figure [5] The negative valued 
quantity is difficult to interpret in terms of the density of 
realistic objects. Thus, the density of the thermal (real) 
vortices needs to be defined using a nonstandard way 
like it is done for the particlelike objects (|4"5| in Ref. [g|. 
Note, that in our simulations we do not discriminate the 
thermal (wrapped) and the UV vortex components. 




T/T 



FIG. 5: (Color online) The difference between the vortex den- 
sity at finite temperature and at zero-temperature (physical 
units) . 

One can observe a few interesting features of the dif- 
ference in the vortex densities shown in Figure [5] 

• The difference in the vortex densities is positive in 
the confinement phase, T < T c . 

• This quantity is negative in the deconfinement 
phase up to temperatures of the order of, approxi- 
mately, three critical temperatures, T c < T < 3T C . 
At higher temperatures the quantity will again take 
positive values. 

• The minimum of the thermal density is reached at 
T w 2T C . Note that this is the temperature at 
which the contribution of the magnetic gluons van- 
ishes, see Eq. (fl?5)) and Figured] 



2. Comment on monopoles vs vortices 

The fact that the minimum of the vortex density is 
reached at T rs 2T c is intriguing. As we have dis- 
cussed above, magnetic vortices are related to magnetic 
monopoles. According to Ref. Q the condensed (in the 
confinement phase) state of the monopoles melts into a 
monopole liquid at T = T c which is evaporated into a 
monopole gas at T « 2T C , at which the vortex density 
reaches its minimum. The coincidence of the temper- 
atures may be not accidental. An illustration of the 
monopole states in the finite-temperature Yang-Mills 
theory is provided in Figure [SJ 



condensate 



liquid 



gas 
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FIG. 6: (Color online) The states of the monopole matter in 
Yang-Mills theory at finite temperature according to Ref. [g|. 



The existence of the monopole liquid close to the crit- 
ical point in the deconfinement phase was confirmed by 
accurate lattice data of Ref. [34[ according to an analysis 
presented in Ref. [35| • The appearance of the monopole 
gas at even higher temperatures - discussed earlier in an- 
alytical terms in Ref. [8| - was demonstrated numerically 



17] 



in Ref. [36J . We refer the interested reader to Ref. 
for a short review of this subject. 

At very high temperatures the general behavior of the 
vortex density can be predicted from dimensional reduc- 
tion arguments. According to the dimensional reduc- 
tion, at high temperatures nonperturbative physics is 
controlled by three-dimensional magnetodynamics which 
corresponds to zero Matsubara frequency of the origi- 
nal four-dimensional theory. The vortex trajectories be- 
come static, and the vortex density becomes equal to the 
density of the vortex lines in the 3d space. The mag- 
netodynamical quantities are expressed in terms of the 
corresponding powers of the gauge coupling in the di- 
mensionally reduced theory, 



9UT)=Tgl d {T), 



(49) 



where g^d is the running coupling of the original four- 
dimensional theory calculated at the scale T. 

Thus, density of vortices in the high-T limit should be 
described by the formula, 



p(T) = C vort g\ d (T) ex 



T 



logT/A 



T » T c , (50) 



where C vor t is a temperature-independent parameter and 
A ~ Aqcd is a dimensional parameter. 

The power, with which the temperature enters 
Eq. (|50[) . is clear from the dimensional arguments, 
p ~ 0(T 2 ). However, the logarithmic behavior, 
log~" vort (T/A), with avort = 2 is governed by the per- 
turbation theory. It is interesting to note that the same 
power of the logarithm, a « 2, was recently found for the 
Abelian monopoles in the high-temperature phase of the 
Yang-Mills theory [34|. This fact is amusing, because for 
the monopoles, which are particlelike objects, the natural 
power of the logarithm would be a mon = 3, Ref. [6, 8]. 
The observed difference allowed the authors of Ref. [34[ 
to conclude that the monopoles do not form a free par- 
ticle gas. We support this opinion by noticing, that the 
monopoles are to be related to the magnetic vortices, and 
the vortex dynamics make the monopole properties less 
trivial. This could explain why the observed power of the 



monopole logarithm, coincides with the vortex one, 

amon ~ "vort = 2 . (51) 

3. Expectation value of SU(2) plaquettes 

We have calculated numerically the expectation values 
of the SU{2) plaquettes ^TrC/p at the original config- 
urations and compared them with the vortex-removed 
configurations. The comparison was done both at zero 
and finite temperature. The plaquette expectation values 
are crucial for the determination of contribution of vor- 
tices into the expectation value of the gluon action, Eqs. 
(gU, flU]), and Jig). Notice, that according to Eq. (fTTij) , 
the larger the plaquette expectation value the smaller the 
action density (and vice versa). 
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FIG. 7: (Color online) The expectation values of the plaquette 
T^Tr Up at zero and at finite temperatures for usual and for 
vortex-removed configurations. 

The expectation value of the SU(2) plaquettes is pre- 
sented in Figure [7] One can immediately observe the 
following interesting features: 

• The plaquette expectation values for configuration with 
removed vortices is much smaller compared to the origi- 
nal unmodified configurations. This feature is valid both 
at zero and at finite temperatures in all phases. 

- Thus, magnetic vortices carry positively valued en- 
ergy density. This fact was already established for 
T = Yang-Mills theory in Ref. (37|. 

• The temperature fluctuations increase the plaquette 
expectation value both at the original configurations and 
at the modified gauge field configurations. 

- Thus, the effect of the temperature is to decrease 
the energy density both at the vortex worldsheets 
and at the space outside of vortices. 

• The influence of the temperature on the plaquette ex- 
pectation values of the vortex-removed configurations is 
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much stronger compared with the observed influence of 
the temperature on the original gauge field configura- 
tions. 

- Thus, the energy densities at the vortex world- 
sheets and outside of vortices receive very strong 
contributions due to thermal configurations. These 
contributions are canceled in the original configura- 
tions while they cannot be canceled in the modified 
configurations due to the removal of the vortices. 

This observation stresses the significance of the vortices 
for the thermodynamics of the gauge system. 

Using the numerical data for the vortex density, for the 
plaquette expectation values both at the original config- 
urations and at the vortex-removed configurations one 
can calculate the contribution of the vortices to the vac- 
uum expectation value of the gluon action density. The 
contribution of the vortices - calculated with the help 
of Eq. (j4Tj) - is presented in Figure M It is clearly seen 
that the effect of the temperature is to decrease the total 
action carried by the vortices. 




FIG. 8: (Color online) The expectation value of the action 
density calculated at the vortex worldsheets (|4ip at zero tem- 
perature (18 4 lattice, the solid line) and at finite temperature 
(18 3 x 4 lattice, the dotted line) vs the gauge coupling /3. 



It is very important to stress, that we have observed 
two different temperature effects which are realized just 
above the critical temperature in the deconfinement phase 
(T C <T<2T C ): 

• According to Figure [H] the higher the temperature 
the lower the contribution of the vortices to the vev 
of the action density. 

• According to Figure [4] the higher the temperature 
the lower the vortex density. 

The knowledge of these facts alone does not allow us to 
conclude what is the primary reason of the dinrinishing of 
the vortex contribution into the total action density. For 
example, the energy density per unit area of the vortex 



worldsheet could be insensitive to the temperature and 
the overall effect could be due to the lowering of the vor- 
tex density with increase of temperature. Other options 
are: the energy density at the worldsheet may moder- 
ately increase (or, decrease) with temperature weaken- 
ing (enhancing) the effect of the dropping vortex density. 
The same question may be addressed at higher temper- 
atures, T > 2T C , at which both the vortex density and 
the thermal plaquette values are increasing functions of 
temperature. 

We study numerically the action density per an ele- 
mentary plaquette belonging to the vortex worldsheet 1 : 



(s P y 



(S t 



_ {T,pPpSp) 

(EpPp) 



(52) 



the action density per an elementary plaquette outside 
of the vortex worldsheet 



sp 



(SpY 



1-P 



<Zp(i-pp)Sp) 

(Ep(I-Pf)) 



(53) 



For the sake of brevity we call the quantity (|52j) as "the 
specific vortex action" because it represents the gluonic 
action density normalized by the area of the vortex world- 
sheet. Analogously, we refer to the quantity (|52"|) as to 
"the specific no- vortex action" . 

The results - which are represented in Figure [9] - are 
quite interesting: 
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FIG. 9: (Color online) The same as in Figure [8] but for the 
action density per an elementary plaquette belonging to the 
vortex worldsheet (|52|l and per an elementary plaquette lying 
outside of vortex worldsheets (15311 . 



1 For the sake of completeness one may suggest to perform even 
more refined measurement by separating the contributions from 
the ultraviolet and the infrared vortex clusters following Ref. [37|1 . 
In our approach we have not done this refinement treating both 
the UV and IR vortices on the same footing. 
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• In the low temperature region of the deconfinement 
phase (T c < T < 3T C ) the specific vortex action is 
larger compared with its value at T = 0. 

• In all phases the specific vortex action is larger com- 
pared with the total action and, consequently, to 
the specific no-vortex action. 

• The specific vortex action is very sensitive to the 
temperature variations compared with other action 
densities. 

The last observation again indicates that the vortices 
play a very important role in the thermodynamics of 
Yang-Mills theory. We discuss this issue in the next Sub- 
section. 



C. Center vortices and trace anomaly: 
numerical results 

The lattice calculation of the trace anomaly requires a 
careful renormalization. The renormalization of the total 
gluonic anomaly is quite straightforward ()20|1 : in order to 
cancel zero-point oscillations the expectation value of the 
517(2) plaquette calculated at zero-temperature should 
be subtracted from the finite-temperature plaquette ac- 
tion. According to a common wisdom, both these quan- 
tities contain the same UV divergencies, and the subtrac- 
tion leaves us with the UV-finite thermal contribution 2 . 

The regularization of the vortex contribution is less 
straightforward. According to Eq. (|40|) the gluonic action 
can be split into two pieces corresponding to the vortex 
contribution, (S P ) VOTt and to the rest, (5 P ) novort , both at 
zero and at finite temperature. However, the quantum 
fluctuations of the gluonic fields at the vortex and outside 
of vortex should not, in general, be the same because the 
vortices are infinitely thin objects possessing very strong 
gluonic fields [37]. Taking into account that the vortex 
density is a lively function of temperature, Figure [5j we 
conclude that the UV divergencies of, say, the action den- 
sity <5 P ) vort at T = and T > should not, in general, 
match. Therefore, the renormalized finite-temperature 
value of ( l S , p) novort by a subtraction of the same quantity 
at T = should give, in general, a UV-divergent result. 

In order to illustrate a possible failure of a naive renor- 
malization, we take the lattice expression for the trace 
anomaly (|20[) . and separate the vortex contribution from 
the rest as follows: 



(Sp) T -(S P ) = [(S P ) 



voit 
T 



{(Sp)£ ovort - [1 



-p(T)(Sp)o] 
p(T)](Sp)o} 



(54) 



Note, however, the cautionary remark in Section IIVB II about 
the renormalization of the vortex density: the naively "renor- 
malized" vortex density turns out to be negative, Figure [5] 



Here we used Eq. ([40]) which implies that the vev of 
the action can exactly be split into the vortex contri- 
bution, (5p) vort , and the contribution outside of vor- 
tices, (Sp) novort . This renormalization subtracts from 
the vortex-originated (5p)™ rt the appropriate amount 
of the T = gauge action, (Sp)q. The coefficient of pro- 
portionality is given by the vortex density p: each vortex 
plaquette at T > is regularized by the plaquette action 
at T = 0. The same type of the regularization is done 
for the no- vortex contribution, {Sp} nomrt . 

According to Eqs. ([2"0]) and ([54"]) the trace anomaly can 
then be split into the two regularized contributions 



9{T) 



gvort 
naive 



(T) 



gnovort 
naive 



en. 



(55) 



where 



nvort 
naive 



(T) 
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(56) 
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naive 



(T) 
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6 if 



d(3(a) 
d log a 



i I C \ novort 

{{bp) T 



[l-p(T)](Sp) }.(57) 



The results for the naively regularized quantities 
and (|57p are presented in Figure [TU] (yet another simple 
regularization was implemented in our earlier investiga- 
tion in Ref. [23]). We see from FigureflOlthat the vortices 
provide a large negative contribution while the space out- 
side of vortices is characterized by a large positive- valued 
contribution. The sum of the two gives the total trace 
anomaly. The naively regularized vortex and no-vortex 
results are not sensitive to the phase transition and they 
continue to rise as the temperature is lowered. Thus, we 
conclude the regularization schemes (J56I) and ([57]) are not 
appropriate. 

It seems that the right quantity which characterizes the 
thermal contribution of the vortices to the action density 
is the specific vortex action ([52]) . i.e. the vortex action 
per unit area of the vortex worldsheet. The regulariza- 
tion (i.e., the subtraction of the zero-point fluctuations 
at zero-temperature) should be done accordingly: the ac- 
tion per vortex plaquette at T > should be regularized 
by the action per vortex plaquette at T = 0: 



(sp)rm = ( s p)t - (sp)o 



(58) 



where the specific vortex action (sp) is defined in 
Eq. (52). 

One can define the specific contribution $ vort of the 
vortices into the trace anomaly, or, the contribution of 
the magnetic vortices to the trace anomaly counted per 
unit area of the vortex world surface at given tempera- 
ture. Formally, the definition of this quantity is 



$vort _ 



de v 



d Area 



(59) 



where "Area" denotes the area of the vortex worldsheet 
at given temperature. The definition of this quantity can 
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FIG. 10: (Color online) Naively regularized contribution of 
vortices into the trace anomaly (|56|) vs temperature (the data 
points are connected by the dotted line). The contribution to 
the anomaly coming from the space outside of vortices (|57[1 , 
and the total gluonic trace anomaly (|20[) are also shown (the 
dashed line and the solid line, respectively). 



be written using Eqs. (j2"0)) and ([5 
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= QLt 



80(a) 
d log a 



((sp)t - (sp)o) • (60) 



-TO 




anomaly from vortices, t 

per vortex area (lattice units) 



T/T 



FIG. 11: (Color online) The specific anomaly, Eqs. (|60|l 
(|59[) . coming from the magnetic vortices (lattice units). 



and 



We show the specific contribution of the magnetic vor- 
tices to the trace anomaly in Figure [Til There are a few 
interesting features of this quantity: in the deconfinement 
phase T c < T < 3.5T C the specific anomaly is negative. 
In the deeper deconfinement phase as well as in the con- 
finement phase the specific trace anomaly is positive. At 
the phase transition the specific anomaly is close to zero. 
All these features are in qualitative agreement with our 
preliminary calculations reported in Ref. 23]. 

The negative (ghostlike) value of the vortex-originated 
anomaly is not unexpected according to the theoretical 



analysis of Ref. [38] . The vortices can be treated as non- 
abelian strings with a nontrivial worldsheet dynamics as- 
sociated, in particular, with the monopoles, localized on 
worldsheets of the magnetic vortices. 

In physical units the specific vortex contribution to the 
anomaly is represented in Figure 1121 The vortex contri- 
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FIG. 12: (Color online) The same as in Fig. llll but in physical 
units. 

bution is numerically very large. This fact agrees quali- 
tatively with a singular nature of the distribution of the 
action density around vortex worldsheets [371 ] . 

Finally, the total regularized contribution of the vor- 
tices into the trace anomaly, 



mrt (T)=p{T)$ vort (T) 



(61) 



is represented in Figure [131 The scale of the vortex- 
originated contribution is useful to compare with the to- 
tal gluonic anomaly, Figure [2] Some features of the total 
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FIG. 13: (Color online) The contribution of the magnetic 
vortices into the trace anomaly, given by Eqs. (|6ip and (|60p . 



vortex contribution to the anomaly (Figure I13[) are very 
similar to the specific trace anomaly coming from vor- 
tices (Figure [T2"|) : in the confinement phase, T < T c , the 
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vortex trace anomaly is a positive quantity while in the 
deconfinement phase this quantity becomes negative in 
the region of temperatures T c < T < 3T C . As for higher 
temperatures, T > 3T C , our data does not allow us to 
discriminate between the two possibilities: the contribu- 
tion to the trace anomaly coming from the vortices may 
be either positive or negative, and the ratio 6* vort /T 4 in 
this region may be very small. Note that due to the def- 
inition of the vortex-originated anomaly, Eqs. ([61"]) and 
([rJTJ]) . this quantity should be equal to zero at T = 0. At 
the phase transition the anomaly is close to zero too. 

It is appropriate to notice that at finite temperature 
some vortex trajectories may wrap with respect to the 
temporal (temperature) direction while other vortices 
may have zero wrapping index. Following Ref. [6[ we 
suggest here that the wrapped vortices have a direct re- 
lation to the thermodynamics of the system while the un- 
wrapped vortices bear features of the zero-temperature 
theory only. Thus we expect that the vortices with dif- 
ferent wrappings with respect to the temporal direction 
may give different contributions to the gluonic anomaly. 
In our simulations we do not discriminate between the 
vortices with different wrapping properties. 

Certain qualitative features of the vortex anomaly, Fig- 
ure 113] may in principle be related to certain geometri- 
cal (orientational) properties of the vortex worldsheets. 
Indeed, it is well known [lg, [2JJ, that in the deep con- 
finement phase, T -C T c , the worldsheets of the mag- 
netic vortices have no preferred orientations in the lat- 
tice spacetime because there is practically no difference 
between spatial and temporal (temperature) directions. 
However, as the temperature increases the preferences 
do appear: the asymmetry ratio between the densities 
of the spacelike (*-Py) and timelike (*Pu, hj — 1,2,3) 
vortex worldsheets slowly rises in the confinement phase 
up to a critical point [13 . l2l| . In the deconfinement 
phase the ratio falls down drastically since at T > T c 
the vortices are getting static according to the dimen- 
sional reduction arguments. On the other hand, we know 
that the chromomagnetic (spacelike, 9m ~ Trf7p i3 ) part 
of the anomaly is always smaller than the chromoelec- 
tric (timelike, 9e ~ Tr Up u ) part, Figured] Since the 
spacelike vortices (*P%j) pierce timelike {Pa) plaquettes 
and vise versa, the slow rise of the spacelike portion of 
the vortex worldsheets in the confinement phase leads to 
an increase of the role of the timelike plaquettes in the 
vortex-dominated anomaly. 

The timelike (chromoelectric) part of the bulk anomaly 
is bigger than the chromomagnetic part according to 
Figure [2] If the same is true for the anomaly at the 
vortex worldsheet then the gradual orientation of the 
magnetic vortices towards spatial directions with the 
rise of temperature should lead to an increase of the 
vortex-originated anomaly in accordance with our ob- 
servations (Figure [13]). As the system passes the de- 
confining point, the timelike vortices become more fa- 
vored compared with the spacelike ones and therefore 
the spacelike chromomagnetic plaquettes become domi- 



nant in the vortex-originated contribution. Thus, taking 
into account the relation 6m < 9e and using the geo- 
metrical arguments, one can qualitatively explain a drop 
of the vortex-originated anomaly just above the critical 
temperature, Figure [T3l 

It is also important to stress that just above the phase 
transition the trace anomaly coming from the vortices 
(Figure [T3"]) is approximately of the same order as the 
total gluon anomaly (Figure [2]). This fact highlights the 
significance of the vortex degrees of freedom to the ther- 
modynamical behavior of the gluon plasma. 

Our numerical calculations do not allow us to trace 
the behavior of the vortex-originated trace anomaly in 
the confinement phase at low enough temperatures, T < 
0.8T C . As this quantity is still not zero at low tempera- 
tures, we cannot calculate the contribution of the vortices 
into the pressure ([12]) and into the energy density (fl3]) 
because they involve integration over the low tempera- 
ture region. More refined and accurate calculations are 
needed to address this issue properly. 

A qualitative extrapolation of our data into the low- 
temperature domain of the phase diagram allows us to 
suggest the behavior of the vortex contribution to the 
trace anomaly as it is visualized in Figure [T¥l 



0.2 
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FIG. 14: (Color online) An expected behavior of the contribu- 
tion of the magnetic vortices into the gluonic trace anomaly 
(the solid line). The dotted line represents another high- 
temperature behavior that cannot be excluded by the avail- 
able data. 



V. SUMMARY OF NUMERICAL RESULTS 

We observe the intriguing relevance of magnetic vor- 
tices to the thermodynamics of Yang-Mills theory. A 
brief summary of the vortex properties is as follows: 



Vortex density, Figure [5] 

• T < T c . In the cold region of the confinement phase 
p(T) is an increasing function of temperature T. 
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As the temperature gets higher, the increase of the 
density turns into a sharp fall. 

• T = T c . At the phase transition the vortex density 
is close to its zero temperature value, p(T c ) ~ p(0). 

• T c < T < 2T C . Right above the critical temper- 
ature the vortex density continues to fall and the 
minimum of the vortex density is reached approx- 
imately around T rs 2T c . At this temperature the 
magnetic contribution to the trace anomaly van- 
ishes, Figure [2] 

• 2T C < T < 3T C . The vortex density start to 
rise, and at T « 3T C the density reaches its zero- 
temperature value again, p(3T c ) rs p(0). 

• T > 3T C . In this domain of the temperatures the 
density is a positively defined rising quantity. 

• T 3> T c . At very high temperatures the vortex 
density should follow the law of Eq. ([50]) . 

Our results on the vortex properties are in a qualitative 
agreement with other studies of the vortex properties at 

finite temperature [H E3, HH, H^] • Thc risc and fal1 of 
the vortex density with the increase of temperature in the 
confinement phase should be understood as a result of the 
competition of two opposite processes: the thermal fluc- 
tuations (i) excite the vortex media leading to creation of 
additional virtual vortex loops from the vacuum; (ii) de- 
stroy the chaotic percolation of the vortex cluster in the 
temporal direction making the vortex trajectories more 
static. The percolation of the vortex trajectories in the 
temporal direction is not a well defined quantity because 
of the finiteness of the extent of the Euclidean spacetime 
in this direction. Nevertheless, at low temperatures not 
all the vortex loops wind around this direction and the 
increase of the temperature tries to break the loops into 
pieces. 



Vortex contribution into the trace anomaly, 
Figures [13] and [14] 

We found that vortices do contribute to the equation 
of state of the gluon plasma. The vortex contribution 
is illustrated in Figures [13] and [M] Qualitatively, the 
temperature behavior of the vortex trace anomaly is very 
similar to the behavior of the vortex density described 
above. We observe that: 

• in the confinement phase the vortex trace anomaly 
is positively defined, 

• at the critical temperature, T — T c , the vortex- 
originated trace anomaly is consistent with zero, 

• in the low-temperature domain of the deconfinc- 
ment phase (T c < T < 3.5T C ) the trace anomaly 
coming from the vortices is negative, 



• at the high-temperature region of the deconfinc- 
ment phase (T > 3.5T C ) the contribution of the 
vortices to the gluon trace anomaly may either be- 
come positive or stay negative. Anyway, the ratio 
^vort (r)/T 4 is expected to be small in this region, 

• qualitatively, the features of the specific vortex 
trace anomaly (Figures QT] and [T2"]) and the "bulk" 
vortex trace anomaly (Figure I13p are very similar 
to each other, 

• the gluon action density at the vortex trajectories 
(Figure [5]) is much higher compared with the gluon 
action density outside of vortex worldsheets. This 
result is an agreement with previous studies re- 
ported in Ref. [37|, and 

• the value of the gluon action density at the vortices 
is much more sensitive to the temperature varia- 
tions compared with thc gluon action density out- 
side of vortices (Figure [5]). 

One can characterize the very strong contribution of 
the vortices into the trace anomaly and, consequently, to 
the equation of the state of the Yang-Mills plasma as fol- 
lows: on our lattices the vortices occupy just 2-4% of the 
lattice spacetime (Figure 2]) in the deconfinement phase 
while their regularized contribution to the trace anomaly 
(Figure [T3"]) is of the order of the total contribution of all 
gluons in the bulk of the system (Figure [2]) . 



VI. CONCLUSIONS 

The basic idea of this paper is to calculate the con- 
tribution of vortices to the vev of the Euclidean action 
density of the gluonic fields. The action density, which 
is basically the gluon condensate, is related to the ther- 
modynamic quantities of the Yang-Mills theory - such 
as pressure and energy density - via the anomaly of the 
trace of the energy-momentum tensor. A short summary 
of the results is presented in the previous Section. 

We observed that magnetic vortices play an essential 
role in the thermodynamics of Yang-Mills theory. Our 
main result is presented in Figure [T3"l In the temperature 
range T c < T < 3T C the contribution of the magnetic 
vortex strings to the trace anomaly is numerically large 
and negative in sign compared to the total energy and 
the trace of the energy-momentum tensor of the whole of 
the plasma: 

flvort (T) < , (Tr G 2 (T) } vort > , T c < T < 3T C . 

(62) 
Taken at face value, Eqs. (|6"2"]) amount to the observation 
of a ghostlike component of the gluonic plasma. 

The existence of magnetic vortices may strongly influ- 
ence the physics of quarks in the quark-gluon plasma. 
The random chromomagnetic fluxes of the string net- 
works in the gluon plasma may trap, scatter and drag 
the quarks. The chromoelectric and chromomagnetic flux 
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tubes are suggested to arise in the early Glasma regime 
of the evolution of color gauge fields in high energy heavy 
ion collisions [33]. We leave to a future publication [40| 
discussions of the magnetic string effects in thermally ex- 
cited gluon media. 

There are, however, reservations to be made. First, the 
theoretical interpretation of data referring to a compo- 
nent of the plasma is subject to an uncertainty because 
the magnetic component obviously interacts strongly 
with the rest of the plasma. And it is rather the to- 
tal energy that has a direct physical meaning than the 
distribution of it. Thus, a detailed interpretation of the 
data (|62p asks for further efforts on the continuum-theory 
side. Second, it might be worth emphasizing that, obser- 
vationally, the negative sign is due to a decrease of the 
density of the strings (compared with the T < T c case). 
At this moment, it is not clear whether this mechanism, 
behind the signs in Eq. (|62[) is important for the interpre- 
tation. 

In any case, the numerical calculations reveal that 
magnetic strings, which at our lattices occupy (2-4)% 
of the total volume are crucial for the dynamics of the 
whole plasma (in the limit of the vanishing lattice spac- 
ing a — > this fraction tends to zero). The absolute 
value of their contribution to the energy of the plasma is 
of the order of the total value of the energy in the bulk 
of the system. Thus, perturbative calculations, like |T]), 
agree with the data averaged over the whole plasma while 
nonperturbative dynamics is responsible for a highly non- 
trivial sharing of the same energy density between various 
components of the plasma. A remote analogy is a quark- 



resonance duality. While simple quark graphs describe 
well the cross section averaged over a large energy inter- 
val, resonances are responsible for a more local-in-energy 
structure. It is becoming obvious that a detailed under- 
standing of the plasma properties is not possible without 
clarifying the role of its magnetic component [6j . 
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